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ABSTRACT 

We present near infrared photometry of the Type la supernova 1999ee; also, 
optical and infrared photometry of the Type la SNe 2000bh, 2000ca, and 200 Iba. 
For SNe 1999ee and 2000bh we present the first-ever SN photometry at 1.035 /xm 
(the y-band). We present K-corrections which transform the infrared photom- 
etry in the observer's frame to the supernova rest frame. Using our infrared 
K-corrections and stretch factors derived from optical photometry, we construct 
JHK templates which can be used to determine the apparent magnitudes at 
maximum if one has some data in the window —12 to +10 d with respect to 
T{Bmax)- Following up previous work on the uniformity of V minus IR loci of 
Type la supernovae of mid-range decline rates, we present unreddened loci for 
slow decliners. We also discuss evidence for a continuous change of color at a 
given epoch as a function of decline rate. 

Subject headings: supernovae, photometry; supernovae 

1. Introduction 

Type la supernovae (SNe) are well known to the cognoscenti as standardizable candles 
for determining extragalactic distances. They are considered to be the explosions of CO white 
dwarfs in binary systems. Using optical light curves, PhiUips (1993), Hamuy et al. (1996a), 
and PhiUips et al. (1999) established the relationship between the absolute magnitudes at 
maximum and the decline rates of these SNe. The characteristic parameter is Ami5(i?), the 
number of magnitudes that the SN declines in the first 15 days after S-band maximum. 
Riess, Press, & Kirshner (1996), Riess et al. (1998), and Jha, Riess, & Kirshner (2004) 
have elaborated three versions of the Multicolor Light Curve Shape (MLCS) method. This 



^^The method actually uses the B-, V- and /-band light curves to obtain a characteristic value of /^m.i^{B). 
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method now correlates the light curve shapes in U BVRI with the absolute magnitudes. In 
the MLCS method the characteristic parameter is A, the number of y-band magnitudes 
that a Type la SN at maximum is brighter than (A < 0), or fainter than (A > 0) a fiducial 
object of My(max) = —19.47 on a scale where Ho = 63 km s^^ Mpc^^ . Finally, the stretch 
method of Perlmutter et al. (1997) is used to stretch or shrink B- and V-band light curves 
in the time axis to match a fiducial. Since brighter SNe are more slowly declining, the 
stretch factor is less than 1.0 for the intrinsically brighter SNe, and greater than 1.0 for the 
intrinsically fainter ones. Typically, Amis (5) = 1-1 corresponds to A = and stretch factor 
s = 1.0. 

While the optical light curves of most Type la SNe fit into some kind of an ordering 
scheme, with the marked increase of SN discoveries in the past few years, more and more 
objects are being found which do not fit the patterns. SN 1999ac was a slow riser, fast 
decliner (Labbe et al. 2001, PhiUips et al. 2003). SN 2000cx was a fast riser, slow decliner 
and had unusually blue B — V colors after T(S,„„a:) + 30 d (Li et al. 2001, Candia et 
al. 2003). Recently, Li et al. (2003) describe the even stranger SN 2002cx, which had 
premaximum spectra like SN 1991T (the classical hot Type la SN), a luminosity like SN 
1991bg (the classic sub-luminous object), a slow late-time decline, and unidentified spectral 
lines. These "pathological" examples should motivate breakthroughs in the modeling of 
Type la SNe, allowing us to understand better the more normal objects. 

Following the pioneering papers of Elias et al. (1981, 1985) and the photometry of 
SN 1986G (Frogcl ct al. 1987), essentially no infrared light curves of Type la SNe were 
obtained until the appearance of SN 1998bu. (See Meikle 2000 for a compilation of all 
the IR photometry available three years ago.) In 1999 two of us (MMP and NBS) started 
regular observing campaigns at Las Campanas Observatory and Cerro Tololo Inter- American 
Observatory to obtain optical and infrared light curves of supernovae, primarily Type la's. 
For this we have used mainly the LCO 1-m Swope telescope and the CTIO 1-m Yale-AURA- 
Lisbon-Ohio (YALO) telescope. 

In this paper we report optical and infrared photometry of the Type la SNe 2000bh, 
2000ca, and 2001ba. We present infrared photometry of SN 1999ee. Previous extensive 
UBVRIz optical photometry of SN 1999ee has been pubhshed by Stritzinger et al. (2002). 

A significant fraction of the data reported in this paper was obtained for the Supernova 



^^The stretch method does not work with i?-and /-band hght curves, except around T{Bmax)i because of 
the shoulder in tlie _R-band hght curves and the secondary maximum in the /-band ones. See Fig. 10 below. 
For roughly 90 percent of Type la SNe, the strength of the secondary hump is correlated with Am.i^(B) 
(Krisciunas et al. 2001), hence with the absolute magnitudes at maximum. 



-4- 



Optical and Infrared Survey (SOIRS), a project organized and carried out in 1999 and 2000 
(Hamuy principal investigator). See Hamuy et al. (2001) for further details. SOIRS also 
included a photographic search for supernovae carried out at Cerro El Roblc (Maza 1979). 
Amongst SOIRS discoveries arc three of the four SNe reported here. Wc note that the optical 
and IR photometry and spectroscopy of SN 1999ee (one of the prime SOIRS targets) is the 
largest such dataset ever obtained for a Type la SN. 

2. Observations 

2.1. Photometric Calibration 

Supernova light curves can be derived from images taken on photometric and non- 
photometric nights thanks to the use of digital array detectors. This works providing that 
calibration has been done on some photometric nights, and it assumes that the exposures are 
long enough for any variations of transparency of the sky covered by the digital array to even 
out. In our experience this works better for optical photometry than infrared photometry. 
In the latter the field of view is usually quite small, and one takes many short exposures 
while dithering the position of the telescope. If the sky is non-photometric, only the center 
of the subsequently constructed mosaic is useable. Thus, one hopes to have an appropriately 
bright field star or two close to the SN so that data from all nights can be utilized. 

One of our concerns is whether observing through clouds affects the photometric colors: 
are clouds greyl Scrkowski (1970) demonstrated that U — B, B — V, and V — R colors were 
affected less than 0.01 mag even though he was observing through up to eight magnitudes 
of cloud cover. Walker et al. (1971) and Olscn (1983) reported similar reassuring results. 

As is common, we have tied the optical photometry of our SNe to the U BVRI standards 
of Landolt (1992a). This allows us to calibrate the field stars near the SNe. However, the 
spectra of SNe are unlike spectra of normal stars, so ideally one would then devise a set 
of corrections to place the SN photometry on a particular photometric system, such as 
that of Bessell (1990). Stritzinger et al. (2002) attempted this with their extensive optical 
photometry of SN 1999ee, but in the end decided not to apply any corrections, citing "rather 
disappointing" results. By contrast, Krisciunas et al. (2003) were much more encouraged 
with their similar endeavors with SN 2001el. We used optical data obtained with the CTIO 
1.5-m telescope, the CTIO 0.9-m telescope, and the YALO 1-m telescope. We were able to 
devise corrections for the B- and V^-bands which placed the photometry on the Bessell filter 
system, tightened up the light curves, and eliminated a known source of systematic error 
in the B — V light curves. However, if we had applied our derived corrections for the R- 
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and /-bands, the light curves derived from imagery obtained on different telescopes would 
have actually been pulled further apart, not tightened up. Hence, such corrections were not 
applied to the R and / photometry. 

For the calibration of the photometry presented here we used a combination of aperture 
photometry within the IRAF^^ environment and analysis using point spread function (psf) 

magnitudes obtained with either DAOPHOT (Stetson 1987, 1990) or John Tonry's version of 
VISTA (Tcrndrup, Lauer, & Stover 1984). Field star magnitudes were typically determined 
via aperture photometry, while instrumental magnitudes of the SNe were PSF magnitudes. 
When necessary we performed image subtraction using a script written by Brian Schmidt 
and based on the algorithm of Alard & Lupton (1998). 

Much of the IR photometry presented in this paper was obtained from imagery obtained 
with the LCO 1-m and LCO 2.5-m telescopes. The LCO 1-m infrared camera utilizes a 
Rockwell 256 x 256 NICMOS-3 HgCdTe array, giving a plate scale of 0.60 arcsec pixel"^. 
The LCO 2.5-m data for SN 1999ee were obtained with the CIRSI camera, giving a plate 
scale of 0.20 arcsec pixer^ The LCO 2.5-m data for SNe 2000bh, 2000ca, and 2001ba were 
obtained with a camera containing a similar Rockwell array to that in the 1.0-m camera, 
but with a plate scale of 0.35 arcsec pixel"^. 

A majority of the J- and i/-band of SN 1999ee was obtained with the optical and 
infrared camera ANDICAM on the YALO 1-m telescope. The infrared camera contained a 
1024 X 1024 HgCdTe array by Rockwell, giving a plate scale of 0.22 arcsec pixel"^. The 
optical camera in ANDICAM, which was used for the SN 2001ba observations presented here, 
utilized a Loral 2048 x 2048 CCD, giving a plate scale of 0.30 arcsec pixel"^. 

It was with the LCO 1-m telescope that Persson et al. (1998) established their system 
of IR standards. Thus, no corrections need to be made to our LCO photometry to place it 
on the Persson system. For the reduction of the YALO photometry to the Persson system 
we used color terms derived from observations of Persson standards made in February 2000. 
For YALO J the color term is —0.043 ± 0.005, scaling J — H. For YALO H the color term is 
+0.015 ± 0.005, scaling J-H. For YALO K the color term is -0.003 ± 0.005, scaling J-K. 
For the case of our YALO IR photometry of SN 1999ee, we also used spectra of SN 1999ee 
(Hamuy et al. 2002) and appropriate filter profiles to correct the YALO photometry for the 
non-stellar spectral energy distribution of the SN. 



^''iRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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In the case of z-band photometry, Hamuy (2001, Appendix B) gives synthetic 2;-band 
magnitudes of 20 spectrophotometric standards (Stone & Baldwin 1983), quoting uncer- 
tainties of ih 0.02 mag. For more information on the z-band filter see Schneider, Gunn, & 
Hoessel (1983). We made UBVRIz observations of the SN 2000bh and 2000ca fields from 
29 January to 1 February 2003, along with some of these spectrophotometric standards and 
Landolt (1992a) fields, to calibrate the 2;-band photometry of the SNe. (Additional BVRI 
cahbration was obtained from images of May 2000, when these two SNe were active.) Details 
of the z-band calibration are given in Appendix A. 

The y-band at 1.035 fim is a new photometric band whose prime advantages are de- 
scribed by Hillenbrand et al. (2002). It is a relatively clean atmospheric window. F-band 
photometry should be relatively independent of the specific IR detector, relatively insensi- 
tive to differences in transparency vs. wavelength and nightly variations of water vapor at a 
given site. For a tabulation and graphical representation of the F-band filter transmission, 
see Table 1 and Fig. 1, respectively, of Hillenbrand et al. (2002). 

If we derive F-band values for Vega, Sirius, and the Sun (see Appendix A of Krisciunas 
et al. 2003), we obtain the following hnear transformation: 

{Y - Ks) = -0.013 + 1.614 (J, - Ks) , (1) 

where Jg ("J-short") and Kg ("K-short") are the equivalent magnitudes in the system of 
Persson et al. (1998). Many of the Persson et al. standards have J, — colors in between 
those of Vega and the Sun, so we may rely on interpolation for the presumed Y — Ks colors. 
Of various Persson et al. stars observed by us in the F-band, we note that P9183 has Y = 
12.459, according to Eq. 1, while Hillenbrand et al. (2002) give Y = 12.296 ± 0.030, which is 
0.16 mag (5.4-0") brighter than our interpolated value. For P9155 Eq. 1 predicts Y= 12.285, 
while Hillenbrand et al. (2002) give Y = 12.363 ± 0.05, which is 0.08 mag fainter or l.S-a 
different. 

We derived the F-band magnitudes of the field stars listed in Table 1 using observations 
taken on 6 nights (1999ee) and 3 nights (2000bh), respectively, with 3 to 5 Persson stars per 
night, whose F-band magnitudes we derived using Eq. 1, rather than relying on only one 
Persson star per night whose F magnitude happens to be given by Hillenbrand et al. Our 
F-band photometry of SNe may be subject to systematic errors, on the order of 0.03 mag, 
perhaps larger. In any case, for the first time we present F-band light curves of Type la 
SNe. 

The field star magnitudes given by us in Table 1 were checked against the corresponding 
2MASS values using SIMBAD. Limiting ourselves to field stars brighter than magnitude 15, 
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beyond which 2MASS data have large uncertainties, we find that our Jj-band magnitudes 
are systematically fainter than 2MASS by 0.014 mag, with a 1-a scatter of ± 0.026 mag. 
For H our field star magnitudes are systematically fainter by 0.028 mag, with an = ± 0.046 
mag. For Kg our field star magnitudes are 0.030 fainter than 2MASS, with (Jk = ± 0.059 
mag. Thus, there are no serious systematic differences between our field star photometry 
and 2MASS values. 



2.2. SN 1999ee 

The Type la SN 1999ee was discovered on 1999 October 7.15 UT (JD 2,451,458.65) 
by Wischnjewsky as part of the SOIRS project and reported by Maza et al. (1999). It 
was located 10 arcscc cast and 10 arcscc south of the core of the SA(rs)bc-type galaxy IC 
5179. CCD photometry reported by Stritzinger et al. (2002) was begun the following day; 
these authors give JD 2,451,469.1 ± 0.5 as the time of 5-band maximum. Our infrared 
observations were begun on JD 2,451,461.6, some 7.5 days before T{Bjnax)- From optical 
data only, using analysis similar to that of Phillips et al (1999), Stritzinger et al. (2002) 
derived Ay = 0.94 ± 0.16 for SN 1999ee. Extensive spectroscopy of this SN is discussed by 
Hamuy et al. (2002). 

While observations of SN 1999ee were ongoing, another SN exploded in the same galaxy, 
the Type Ib/c SN 1999ex. It was discovered by Martin et al. (1999) on 1999 November 9.51 
UT. Stritzinger et al. (2002) were able to detect SN 1999ex in prc-discovery images as early 
as 30 October UT. We will present infrared photometry of this SN in another paper. 

In Table 1 we give infrared photometry of some of the field stars near IC 5179. The 
mean values are based on data taken on six photometric nights. We retain the numbering 
scheme of Stritzinger et al. (2002). Star 3 is clearly the best field star for calibration of the 
SN light curves, as it is close enough to the SN to be on the detector at all times that the 
SN was. Star 1 of Stritzinger et al. was saturated on many of our LCO 1-m frames, so we 
have not used it as a secondary standard. Star 9 is quite faint and star 6 is at the very edge 
of the IR mosaics, in which case it is only useful on photometric nights. Star 17 exhibited 
evidence of some level of variability. Because we decided to tie the photometry of SN 1999ee 
to star 3 only, we have added in quadrature 0.03, 0.02, 0.02, and 0.02 mag, respectively, to 
uncertainties of the resulting Y, Jg, H, and Ks photometry. 

In Table 2 we give IR photometry of SN 1999ee from LCO, and in Table 3 we give 
the YALO data along with corrections to the photometric system of Persson et al. (1998). 
For the LCO and YALO imagery we derived the photometry from PSF magnitudes, using 



-8- 



photometric templates obtained on 2000 November 15-17 UT, more than 394 days after 
the time of 5-band maximum. We assume SN 1999ee had faded away sufficiently after 13 
months that there are no serious systematic errors resulting from the image subtraction. 
Given that SN 1999cc had faded by 4.9 mag in the Y-band by 2000 April 28 and that it 
was not detectable in J^, H, or in late April of 2000, the mid- November 2000 templates 
should be acceptable. 

The YALO J and H photometry has been corrected to the system of Jg and H mag- 
nitudes of Persson et al. (1998) using the method described by Stritzingcr et al. (2002) and 
Krisciunas et al. (2003). We used actual spectra of SN 1999cc described by Hamuy et al. 
(2002). Fig. 6 of Krisciunas et al. (2003) shows the derived filter corrections, which are 
added to the YALO IR photometry to place it on the system of Persson et al. (1998). These 
corrections amount to as much as 0.13 mag, and from 10 to 30 d after T{Bjnax) make the 
YALO J and H magnitudes brighter. 

The derivation of the filter corrections involves the construction of effective filter trans- 
mission profiles, each of which includes (as a function of wavelength) the atmospheric trans- 
mission, the nominal filter transmission, aluminum refiections for the telescope and within 
the instrument, window transmission for the instrument, dichroic transmission (if the in- 
strument contains one), and the quantum efficiency of the detector. The effective filter 
transmission is meant to mimic the natural system of the observations. This natural system 
is transformed to the standard system of Persson et al. (1998) using color terms measured 
elsewhere. 

As a consistency check, the J-band color term from synthetic photometry of Vega, 
Sirius, and the Sun is —0.037 (scaling J — H), which compares well to the actual color term 
from observations of Persson et al. (1998) standards of —0.043. For the H-hand the color 
term from synthetic photometry is +0.020, also scaling J — H, which compares well with 
the actual value from observations of stars of +0.015. These tests show that our model 
bandpasses provide a reasonable match to the actual YALO bandpasses. 

In Fig. 1 we show near-infrared light curves of SN 1999ee. In all infrared bands the 
time of maximum was about 4 days prior to T{Bmax)- 

The reader will note that the corrected YALO photometry still does not quite agree with 
the LCO photometry. From about 15 to 20 days after T{Bmax) the corrected YALO values 
arc still fainter than the LCO values. We discovered that we could reconcile the two sets of 
data if we shifted the YALO J filter 225 A to longer wavelengths and shifted the YALO H 
filter 200 A to shorter wavelengths. The synthetic color term for H (+0.035) would still be 
in reasonable agreement with the actual value (+0.015), given uncertainties of 0.01 to 0.02, 
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but for the J-band there would be considerable disagreement (—0.112 vs. —0.043). This 
shows that a simple wavelength shift does not account for the differences in the photometry 
and that the shape of the model bandpasses must be different than the actual ones, or that 
the spectrophotometry of SN 1999ee is in error. The conservative user of our photometry 
should give greater weight to the LCO photometry, since it was primarily taken with the 
very telescope used to set up the photometric system of Persson et al. (1998). 

2.3. SN 2000bh 

SN 2000bh was discovered on 2000 April 5.23 UT (JD 2,451,639.73) by Antezana as 
part of the SOIRS project. The discoverey was reported by Maza et al. (2000a). SN 2000bh 
was located 8 arcsec west and 11 arcsec south of the nucleus of ESO 573-14. A spectrum 
taken on April 7 UT by L. Ho revealed it to be a Type la SN near maximum light. 

Fig. 2 shows the field of SN 2000bh and the nearby field stars. 

Infrared photometry of some of the fields stars near SN 2000bh is to be found in Table 
1. Optical photometry of the nearby field stars is given in Table 4. Optical and infrared 
photometry of SN 2000bh is given in Tables 5 and 6, respectively. Because stars 2 and 3 
were at the edges of the IR mosaics, we have tied the IR photometry to stars 1 and 11 only. 

In Fig. 3 we show the optical and IR hght curves of SN 2000bh. We have added BVI 
fits derived using the Amis method (Phillips et al. 1999). 

We find from analysis similar to Phillips et al. (1999) that Ami5(5) = 1.16 ± 0.10, 
T{Bmax) = JD 2,451,635.29 ± 0.28, E{B - V)host = 0.02 ± 0.03, and m-M = 35.01 ± 0.08 
on the Freedman et al. (2001) distance scale. 

In Fig. 3 we also show near infrared light curves of SN 2000bh. If the behavior of 
SN 2000bh were hke other Type la SNe, its IR maxima occurred about JD 2,451,632, some 
8 d before the start of our observations. We note the extremely deep dip in the Js-band 
light curve. We also note that the i^-band secondary maximum must have been comparable 
in brightness to the first H-hand maximum; this is similar to the behavior of SN 2001el 
(Krisciunas et al. 2003). 

2.4. SN 2000ca 

SN 2000ca was discovered as part of the SOIRS project by Antezana on 2000 April 
28.18 UT (JD 2,451,662.68) roughly 1 arcsec east and 5 arcsec north of the nucleus of ESO 
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383-32. The discovery was reported by Maza et al. (2000b). A spectrum taken by Aldering 
& Conley (2000) on April 29.3 UT revealed SN 2000ca to be a Type la SN near maximum 
light. From the ratio of Si II lines at 580 and 610 nm, Aldering & Conley suggested that 
this SN was slightly hotter and more luminous than normal; also, that the equivalent width 
(0.02 nm) due to Na D at the redshift of the host galaxy indicated modest extinction by the 
host galaxy. 

Fig. 4 shows the field of SN 2000ca and the nearby field stars. 

Infrared photometry of some of the fields stars near SN 2000ca is given in Table 1. 
Optical photometry of the nearby field stars is given in Table 4. Optical and infrared 
photometry of SN 2000ca is to be found in Tables 7 and 8, respectively. For the most 
part we found that photometry based on image subtraction was very close (0.01 mag) to 
photometry based on PSF magnitudes without image subtraction, but to be on the safe side 
we preferred to derive our results from image subtractions. 

The optical photometry was derived using image subtraction templates obtained with 
the CTIO 1.5-m telescope on 2003 February 1 UT. Infrared photometry based on images 
taken with the LCO 1-m telescope was derived using subtraction templates obtained on 2003 
March 9 UT. Infrared photometry based on images taken with the LCO 2.5-m telescope was 
obtained from PSF magnitudes without the use of Jg- and if-band subtraction templates. 

In Fig. 5 we show UBVRIz optical and JgH light curves of SN 2000ca. We find 
from analysis similar to Phillips et al. (1999) that Ami^{B) = 0.98 ± 0.05, TiBmax) = JD 
2,451,666.22 ± 0.54, E{B-V)host = 0.00 ± 0.03, and m-M = 34.91 ± 0.08 on the Freedman 
et al. (2001) distance scale. 

2.5. SN 2001ba 

SN 2001ba was discovered by Chassagnc (2001) from images taken on 2001 April 27.8 
and 28.7 UT (Julian Dates 2,452,027.3 and 2,452,028.2). It was located 19 arcsec east and 
22 arcsec south of the nucleus of MCG-05-28-1. A spectrum by Nugent & Wang (2000) 
obtained on April 30 UT revealed it to be a Type la SN near maximum light. 

Fig. 6 shows SN 2001ba and the nearby field stars. Optical and infrared photometry 
of the field stars is given in Tables 4 and 1, respectively. The field stars were calibrated 
from observations with the LCO 1-m telescope. All of the optical data of SN 2001ba itself 
were obtained with the YALO 1-m telescope. From 20 nights of YALO data we determined 
mean color coefficients, using our LCO cahbration of the field stars near SN 2001ba as local 
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"standards" . We then used these mean color terms to determine photometric zero points 
for the transformation of the instrumental PSF magnitudes of the SN to standardized values. 
This resulted in the smoothest possible light curves. 

Optical photometry from the YALO 1-m telescope is given in Table 9. IR photometry 
from the LCO 1-m and 2.5-m telescopes is given in Table 10. The optical and infrared hght 
curves of SN 2001ba are shown in Fig. 7. 

We find from analysis similar to Phillips et al. (1999) that Ami5(5) = 0.97 ± 0.05, 
T^{Bmax) = JD 2,452,034.48 ± 0.68, E{B - V)host = 0.04 ± 0.03, and m-M = 35.55 ± 0.07 
on the Preedman et al. (2001) distance scale. 

We have not devised filter corrections to place the YALO optical photometry of SN 2001ba 
on the system of Bessel (1990). As explained in our paper on SN 2001el (Krisciunas et al. 
2003), the B — V colors in the tail of the color curve will likely be "too red" by ^ 0.1 mag 
as a result. This systematic difference was taken into account when calculating E{B — V)taii 
as part of the Amis analysis. 

3. Discussion 

3.1. Determining mcLximum magnitudes in the infrared 

It is well known that the optical light curves of Type la SNc show various patterns. For 
example, if wc order the light curves of different Type la SNe by the decline rate parameter 
Amis (5), the MLCS luminosity parameter A or the Perlmutter et al. (1997) stretch factor, 
we see that the secondary hump in the /-band decreases in strength and occurs closer to the 
time of S-band maximum as we move from slower to faster decliners. Narrower light curves 
in B and V have weaker 7-band secondary humps (Hamuy et al. 1996b). Por the fastest 
decliners such as SNe 1991bg and 1999by there is no secondary hump at all (Garnavich et 
al. 2001, and references therein). Por the i?-band one typically sees a "shoulder" in the light 
curves of the slow decliners, but no shoulder for the fast decliners. 

Ehas et al. (1985), Meikle (2000), Krisciunas et al. (2000), Verdugo et al. (2002), and 
Phillips et al. (2003) have sought templates for the JHK hght curves, with mixed results. 
Candia ct al. (2003) also give some JHK light curves, ordered by Amis (5). Along with the 
data presented in this paper, we can state that the morphology of the infrared light curves 
does not change with the same obvious pattern as the i?-band and /-band light curves. 
Roughly speaking, the time of the secondary J-band hump moves to earlier phases for faster 
decliners. The strength of the dip in the J-band light curves roughly two weeks after T{Bmax) 
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and the strength of the secondary hump following are not a monotonic function of Ami5(S). 
Also, the flatness of the i?-band light curves 10 days after T{Bjnax) is not a simple function 
of Ami5(S). 

Given the available data, we have not yet devised a method of fitting the JHK light 
curves from t — —12 to +60 days with respect to the S-band maximum. However, we show 
here that there are stretchable templates that appear to fit the infrared light curves within 
~13 days of the times of the infrared maxima. 

First it is necessary to correct the infrared photometry from the observer's frame to the 
SN rest frame. See Appendix B for details on the infrared K-corrections. These corrections 
are given in Table 11 and are to be subtracted from the data to correct them for the Doppler 
shifting of the spectral energy distributions of the SNe with respect to the filter profiles. 

Jha (2002, Fig. 3.8) gives new relationships between Ami5(S) and the stretch parame- 
ters for V and B. The reader will note that the stretch parameter s is defined as the factor by 
which one stretches a template to match a V- or S-band light curve. If we seek a template, 
we seek values of to apply to light curves. 

Using values of Ami5(S) from PhiUips et al. (1999), Krisciunas et al. (2001, 2003), 
Strolger et al. (2002), and this paper, we calculated the inverse stretch factors (s"^) from 
the V- and S-band regressions of Jha (2002). averaging the two for each object. The objects 
in question were SNe 1980N, 1986G, 1998bu, 1999aw, 1999ee, 2000ca, 2001ba, and 2001el. 
All of these objects were observed overlapping the times of IR maxima, or shortly thereafter. 
This made it easy to estimate the maximum infrared magnitudes. 

Our templates use "stretched time" as the time coordinate. We transformed the dates 
of the observations to (Julian Date — T{Bmax)) times s^^ and took out the cosmological 
time dilation by dividing by {1 + z), where z is the redshift. 

In Fig. 8 we show the superimposition of the K-corrected, stretched J-band light curves 
with respect to maximum light. While the dispersion increases significantly after t = +10 
d, the stretched light curves show very uniform behavior overlapping the time of maximum. 

In Fig. 9 we show the JHK templates overlapping the time of the infrared maxima. 
The RMS residuals of the third order polynomial fits are ± 0.062, 0.080, and 0.075 mag, 
respectively, for J, H, and K. In Table 12 we give polynomial coefficients so that the reader 
may use the fits shown in Fig. 9 for the extrapolation to maximum light for other supernovae. 

The reader might naturally wonder: does the stretch method work at all for the /-band? 
In Fig. 10 we show the data for SNe 1998bu, 1999aw, 1999ec, 2000ca, 2001ba, and 2001cl 
treated in the same manner as our JHK data. Clearly, the stretch method does work at the 
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time of the /-band maximum. The RMS scatter in the plot is only ± 0.045 mag. 

We are particularly encouraged by results shown in Fig. 9, because they give us a tool 
with which to estimate the maximum infrared magnitudes of other Type la SNe which have 
many fewer data points in the window of stretched time from —12 to +10 days. In particular, 
Krisciunas, Phillips, & Suntzeff (2003) use the extinction-corrected maximum magnitudes of 
the template objects listed above, along with data of several other Type la SNe to produce 
JHK Hubble diagrams. They further show that there are no obvious decline rate relations 
if we plot the resultant absolute magnitudes of the SNe vs. Amis (5). Thus, Type la SNe 
in the infrared may be standard candles when at maximum brightness. 

3.2. Color curves and extinction 

Until recently, extinction corrections for SNe have primarily been calculated from a 
B — V color excess and an assumed value of Ry, where Ay = Ky x E{B — V). Standard 
Galactic dust gives Ry = 3.08 ± 0.15 (Sneden et al. 1978). But there is no universal value 
of Ky. For example, in the case of SN 1999cl Krisciunas ct al. (2000) found that Ry = 1.8. 
Similarly, Hough ct al. (1987) found that Ry = 2.4 ± 0.13 may be appropriate for the dust 
in NGC 5128, the host of SN 1986G. Now that supernovae are regularly observed at optical 
and infrared wavelengths, we can use a variety of color indices to calculate Ay, perhaps with 
fewer worries about systematic errors. 

In our first paper on optical and IR photometry of Type la SNe (Krisciunas et al. 
2000) we constructed unreddened loci of objects with mid-range decline rates. These loci 
were based on very few data points per object. However, the well sampled photometry of 
SN 2001el and color curves based on a generic Hofiich model (Krisciunas et al. 2003, §3.4) 
showed that our original V — H and V — K unreddened loci contain no serious systematic 
errors. This is to say that our empirical color curves and the synthetic color curves based 
on a SN model are consistent with each other. As stated by Krisciunas ct al. (2003), the 
agreement is better than the accuracy of the model. The bottom line is that we believe the 
unreddened loci derived by Krisciunas et al. (2000) allow us to determine Ay to better than 
0.1 mag for Type la SNe of mid-range decline rates. 

Over the past four years we have obtained enough well sampled infrared light curves 
of Type la SNe that we may now construct the unreddened loci for Type la SNe which are 
slow decliners. Also, we can address the question of the "uniformity" of the color curves for 
the slow decliners. Candia et al. (2003) gave prehminary V — J and V — H loci for the slow 
decliners. What follows supercedes that work. 
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Krisciunas et al. (2000, 2001) previously considered objects of relatively low redshift 
z < 0.01). Since we are now investigating objects up to 2 = 0.038, it is necessary first of all 
to apply K-corrections to the VJHK data. For the V^-band data we interpolated the values 
in Table 4 of Hamuy et al. (1993). For JHK we used the values presented here in Table 11. 

If we adopt the values of Aa / Ay = 0.282, 0.190, and 0.114 for the J-, and K-bands, 

respectively, given by Cardelli, Clayton, & Mathis (1989; Table 3, column 5) and assign an 
uncertainty of 20 percent to each ratio (to account for some range of dust properties in other 
galaxies), it follows that: 

Ay = (1.393 ± 0.110) £;(y- J) . (2) 
Ay = (1-235 ± 0.058) ^(y - i^) . (3) 

Ay = (1.129 ±0.029) £;(y- K) . (4) 

In Fig. 11 we show dereddcncd V minus near infrared colors of 1999aa, 1999aw, 1999ce, 
1999gp, 2000ca, and 2001ba, whose decline rates ranged from Ami5(i3) = 0.81 to 1.00.^^ Our 
adopted values of Ay were derived from the optical light curves using the method of Phillips 
et al. (1999), and our values of V minus IR color excesses were derived using Eqs. 2 to 4. 
The data shown in Fig. 11 are corrected for K-terms, extinction and time dilation. 

In Table 13 we give least-squares regressions to subsets of the data shown in Fig. 11. 
The reader should pay particular attention to the range of time over which we feel the fits 
are valid. 

From —8 ^ t ^ +9.5 d the dcrcddcned V — J colors get monotonically bluer. The 
scatter about the line is only ± 0.067 mag. Thus, from the available data, V — J colors of 
slow decliners may be regarded as "uniform" over this range of time with respect to S-band 
maximum. However, at t — +9.5 d the dispersion oiV — J colors increases significantly. 

In Fig. 11 we show as dashed lines the fourth order fits to the V — H and V — K colors 
of eight SNe of mid-range decline rates considered by Krisciunas et al. (2000). Coefficients 
to generate these polynomials are given in Table 14. These loci are consistent with Peter 
Hofiich's modeling presented in our paper on SN 2001el (Krisciunas et al. 2003, §3.4). One 



^^In the case of SN 1999aa we use the updated value (J. L. Prieto, private communication) of Ami5(B) 
= 0.81 ± 0.04 based on the data of Krisciunas et al. (2000) and Jha (2002). 
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obvious thing to note about the slow dechners is that they have V — H and V — K colors 
which are bluer than those of the mid-range decliners. At the time of S-band maximum the 
slow decliners are 0.243 and 0.230 mag bluer, respectively, in these two color indices. 

In the case olV — H colors, a fourth order fit to the dereddened colors of SNe 1999aw, 
1999ee, 1999gp, 2000ca, and 2001ba exhibits a scatter of ± 0.062 mag with a reduced of 
1.9 prior to t = +8.5 d. However, the scatter from +8.5 < t < +27 d is ± 0.15 mag with 
= 12. Thus, V — H colors of slow decliners may be considered "uniform" within 8 or 9 
days of T^B^ax)-! but not afterwards. 

Our dereddened V — K colors of slow decliners do not exhibit any epoch over which the 
RMS scatter is less than ± 0.1 mag. The V — K colors of SN 1999ee, for example, are roughly 
0.3 mag redder than those of SN 2001ba, and these two objects have quite similar i?-band 
decline rates. While the formal scatter of the fourth order fit to the data of SNe 1999aa, 
1999aw, 1999ee, 1999gp, and 2001ba is ± 0.138 mag, we do not think that the available 
data allow us to assert that the V — K colors of the slow decliners are particularly uniform. 
The reader should notice from Table 11 that the K-corrections are quite large for the LCO 
X^-band and a strong function of the redshift. Some of the spread of the V — K colors shown 
in the bottom panel of Fig. 11 may be due to differences in the X-band spectra of the SNe, 
which is to say that our K-corrcctions, based on SN 1999ee but applied to different objects, 
may add systematic error to the results. Sec Figs. 7 and 8 of Hamuy et al. (2002) for a 
comparison of spectra of different SNe. Some of the spread of the photometry must be due 
to the uncertainties of the reddening coefficients in Eqs. 2 to 4. 

We wish to consider if there is evidence of a continuous change of V minus IR colors 
at some epoch. One can often determine a maximum magnitude more accurately than 
the time of maximum. As a result, we show in Fig. 12 the dereddened colors Vmax — Jmaxi 
Vmax — Hjnaxi ^ud Vmax ~ K max ■^'^ Thc objccts cousidcrcd arc SNe 1980N (Hamuy ct al. 1991), 
1981B (Buta & Turner 1983), 1998bu (Suntzeff et al. 1999, Jha ct al. 1999), Phillips ct al. 
2003), 1999ee (Stritzinger et al. 2002), 2000bh, 2000bk (Krisciunas et al. 2001), 2000ca, 
2001ba, and 2001el (Krisciunas et al. 2003). We also consider SNe 1994D (Richmond et al. 
1995), 1999gp and 2000ce (Krisciunas et al. 2001), which had very few data points in the 
[—12, +10] day window. The IR data can be found in the papers cited, in Meikle (2000), and 
in the tables of this paper. 

Since relatively few SNe are discovered early enough to measure their infrared maxima 



^^The reader will note that because the y— band maximum of a Type la SN occurs on average about five 
to six days after the IR maximum, thc colors plotted in Fig. 12 are unphysical in the sense that they do not 
correspond to observations that can be made at a single moment in time. 
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directly, we plot in Fig. 13 the interpolated V — J and V — H colors 6 days after B-band 
maximum. The regression lines weighted by the errors of the points are: 

{V - J)t=+6 = (-1.837 ± 0.148) + (0.570 ± 0.138) x Ami5(5) . (5) 

{V-H)t=+Q = (-1.836 ± 0.169) + (0.650 ±0.157) x Amis (5) . (6) 

The RMS residual of the F — J fit is ± 0.125 mag, with a reduced value of 1.22, while for 
V — H, cr(RMS) = ± 0.146 mag and xl — 1-32. The scatter corresponds to to an uncertainty 
in Ay of ± 0.18 mag, or ± 0.06 mag in E{B — V). This is comparable to the advertised 
accuracy of determining extinction and reddening with the Amis niethod. 

As time goes on our growing database will allow us to make plots analogous to Figs. 
12 and 13 which will include only those objects with small host extinction and whose light 
curves were well sampled at maximum. These two preliminary figures can be interpreted to 
mean that there is evidence for a continuous change of — J and V — H color as a function 
of decline rate. V — K data may eventually show a similar trend. 

Regarding V minus IR loci for the mid-range decliners and the slow decliners, the 
presently available data can be looked at two ways. It is equally valid to average the data 
of several objects over some range of decline rates, or to assert qualitatively that there is 
a continuous change of color as a function of decline rate at some epoch earlier than t = 
+9 d. The idea of continuous color variation as a function of decline rate is more pleasing 
aesthetically, because one would think that cosmic explosions exhibit a continuous range of 
energies, temperatures, and opacities. 

4. Conclusions 

In this paper we presented well sampled light curves for four Type la SNe. This includes 
the first-ever SN data published at 1.03 fira, though the F-band data lack a firm basis in 
calibration. Given that the maxima in the infrared light curves occur typically 3 to 4 days 
before the time of S-band maximum, it is a considerable challenge to obtain light curves 
that cover the IR maxima. We achieved this in the cases of SNe 1999ee, 2000ca, and 2001ba. 

We have now observed a sufficient number of Type la SNe near the time of maximum 
that we can investigate the question of uniform templates for the infrared. If we stretch the 
light curves in the time domain according to stretch factors based on the B- and y-band 
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relationships of Jha (2002, Fig. 3.8), we obtain reasonably uniform templates from —12 

to +10 days after maximum (in "stretched days"). The RMS uncertainties of the fits are 
± 0.062, 0.080, and 0.075 mag in J, H, and K, respectively. This allows us to determine 
the maximum magnitudes of other Type la SNe as long as we have some data in the —12 to 
+10 d window, and the light curves are not obviously "abnormal" (e.g. like SN 2000cx). 

Primarily on the basis of observations carried out at CTIO and Las Campanas since 
the beginning of 1999, we can now draw the V minus IR color curves for Type la SNe 
which are mid-range decliners and slow decliners. We can utilize the present observational 
data, backed up by modeling for the mid-range decliners (Krisciunas et al. 2003, §3.4), to 
determine the extinction suffered by these objects. This is true even for highly reddened 
objects that occurred in galaxies with different dust properties than we measure for the dust 
in our Galaxy. 

The motivation behind using IR data is that extinction is considerably lower than in 
the optical, and if we know the intrinsic V minus IR colors of Type la SNe, then we can 
determine Ay by slightly scaling the V minus IR color excesses. This should give us more 
accurate extinction corrections compared to deriving B — V color excesses and scaling them 
by an assumed value of Ry. However, while we can enumerate unreddened loci for Type la 
SNe of mid-range and slow decline rates, for the slow decliners the "uniformity" oiV — J 
and V — H color is limited to within nine days of the time of S-band maximum, just about 
what we found for the construction of our JHK templates. 

The available data can also be used as evidence that there is a continuous change of 
V minus IR color at t = +6 d vs. the decline rate. The Vmax — Jmax and Vmax — Hmax colors 
are further evidence that Type la SNe are progressively redder as a function of decline rate. 
Given that we already have evidence for a spectroscopic sequence as function of temperature 
(Nugent et al. 1995), the further calibration of the optical-IR color relationships will be 
relevant to SN modeling. 

As large scale surveys of the sky are planned and implemented, we can and should make 
use of the light curve and color templates presented here to derive the maximum magnitudes 
of Type la SNe and correct the data for the effects of dust along the line of sight. Since 
Type la SNe are the most important distance indicator at redshifts z > 0.01, we may further 
sohdify the foundation of the cosmological distance ladder with these useful cosmic beacons. 

We made extensive use of the NASA/IPAG Extragalactic Database (NED). We also 
made use of the SIMBAD database, operated at GDS, Strasbourg, Prance. We thank the 
Space Telescope Science Institute for the following support: HST GO-07505.02A; HST GO- 
08177.06 (the High-Z Supernova Team survey); HST GO-08641.07A was provided by NASA 
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through a grant from the Space Telescope Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
We thank STScI for the salary support for PC from grants GO-09114 and GO-09421. We 
thank Arlo Landolt for many useful discussions. Our late colleague Robert Schommer was a 
strong supporter of the YALO 1-m telescope, with which we have obtained a large amount of 
supernova data. KK thanks LCO and NOAO for funding part of his postdoctoral position. 

A. Calibration of z-hand Photometry 

Hamuy (2001, Appendix B) gives synthetic 2;-band magnitudes of 20 of the spectropho- 
tometric standards of Stone & Baldwin (1983). Using the VRI photometry of 17 of these 
stars given by Landolt (1992b), we obtain the following regression: 

(R-z) = (-0.0416 ±0.0069) + (0.7627 ± 0.0140) (F -/) , (Al) 

with an RMS residual of ± 0.020 mag. The range oiV — I color of this sample is —0.266 to 
+0.769. 

From 10 to 12 November 2001 and 30 January to 1 February 2003 we imaged a number 
of Landolt (1992a) fields in VRIz using the CTIO 1.5-m telescope and also imaged some of 
the Stone-Baldwin standards. In Fig. 14 we show the R — z vs. V — I colors. 

The solid line in Fig. 14 is that of Eq. Al above. Over a similar color range of the 
Stone-Baldwin standards, the RMS scatter of R — z colors of the Landolt (1992a) standards 
is ± 0.039 mag. Some of this larger scatter is due to three stars whose points he noticeably 
above the fine: Rubin 149C, SA95-96, and Rubin 152E. We offer no explanation for these 
deviations, but note that these three stars have B — V colors like those of stars of spectral 
type A. 

Somewhere around V — I — 1.00 there is a change of slope. From the points with 
V-I> 0.97 we find: 

(R-z) = (-0.0196 ± 0.0289) + (0.7191 ± 0.0217) (V - I) , (A2) 

with an RMS residual of ± 0.030 mag. 

Naturally, if one wishes to calibrate z-band photometry, it is best to use the Stone- 
Baldwin standards, but with the three exceptions noted above, there is a rather tight rela- 
tionship between the R — z and V — I colors of Landolt (1992a) standards. 



-19- 



B. Infrcired K- correct ions 

The K-corrections for the near-infrared bands JgHKs have been calculated in the man- 
ner described in Hamuy et al. (1993). The K-correction is given with the usual sign conven- 
tion: mi{z = 0) = 'mi{z) — Ki{z) where mi{z) is the magnitude as observed at the telescope 
and mi{z — 0) is the corrected magnitude as if it were observed in the frame of the super- 
nova. Optical K-corrections are typically positive and an increasing function in z because 
the optical flux of most sources increases towards the red. The K-correction vanishes for 
Fx oc independent of the filter. Since the near-IR flux tends to decrease as function 
of wavelength, most of the near-IR K-corrections are negative. The K-corrections are not 
as monotonic in z as in the optical because large nebular emission features dominate the 
continuum flux. 

These K-corrections are based on 11 spectra of a single supernova, SN 1999ee (Hamuy 
et al. 2002), spanning days —9 to +41 from T{Bmax)- The spectra were corrected to zero 
heliocentric velocity from the observed velocity of f = 3498 km s^-*^. Each spectrum was 
inspected by eye. For the regions between the J/ H and H/K bands where there is significant 
water absorption, a simple hnear interpolation was made. In some cases, we extended the 
spectrum on a given date by logrithmic interpolation between the spectra bracketing the 
given spectrum in time. This was generally done to extend a spectrum slightly to the red in 
order for the spectrum to overlap "unredshifted" filter curve. 

The filter transmission functions were calculated from the filter curves given in Persson 
et al. (1998) for the J^, H, and Kg filters. These filter functions were muhtplied by a 
standard atmosphere at airmass = 1 (to introduce the saturated telluric absorption), two 
aluminum reflections (to mimic the telescope mirrors) , three aluminum reflections (to mimic 
the internal reflective optics), and a typical quantum efflcieny of a Rockwell HgCdTe detector. 
The remaining optical element — the Dewar window — was ignored. The product of these 
curves defines the natural system of the Persson et al. standards. 

Because these K-corrections are based on a single supernova, these tables should be 
used with caution for supernovae that have significantly different spectral features than 
SN 1999ee. 
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Table 1. Infrared Photometry of Field Stars 
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(0 
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14 
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(0 
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(0 
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16 
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(0 
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13 
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(0 
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4 
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12 
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(0 
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12 


003 


(0 
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(0 
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6 
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^0 
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14 
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(0 
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14 
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(0 
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16 
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^0 
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16 
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(0 
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16 
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(0 


047) 





4 


4 3 




20 




15 


862 


;o 
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15 


460 


(0 


037) 


15 


455 


(0 


025) 





4 


4 3 



^Number of observations that were used to determine the weighted means, for filters F, 
Js, H, and Ks, respectively. 
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Table 2. LCO Infrared Photometry of SN 1999ee 
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V 
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15 
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(0.026) 




2 


478 51 






15 


716 fO 0241 










2 


479.52 










15 


411 
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> 19.00 




> 18.64 


> 17.00 


1 


662.86 


19.64 


(+0.26.-0.21) 














1 



^Julian Date minus 2,451,000. 
•^1 = LCO 1-m. 2 = LCO 2.5-m. 

^The data of Julian Date 2,451,661 are 3-a upper limits. 



-25- 



Table 3. YALO Infrared Photometry of SN 1999ee and Filter Corrections^ 
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60 


16 
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(0 
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-0 
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15 
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(0 
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033 



-26- 



Table 3 — Continued 

517.54 16.887 (0.040) -0.131 16.067 (0.044) 0.033 

519.54 17.129 (0.047) -0.131 16.128 (0.054) 0.033 

521.59 17.263 (0.063) -0.131 16.208 (0.070) 0.033 



^The photometry corrections of columns 3 and 5, when 
added to the data given in columns 2 and 4, respectively, 
place the YALO JH data on the JgH system of Persson et 
al. (1998). These corrections were derived from spectra of 
SN 1999ee (Hamuy et al. 2002) and appropriate filter trans- 
mission functions. 

Julian Date minus 2,451,000. 



Table 4. Optical Photometry of Field Stars 
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Table 4 — Continued 



14 ••• 16.150 (0.024) 15.254 (0.016) 14.761 (0.022) 14.276 (0.023) 

17 ■■• 14.767 (0.023) 13.952 (0.016) 13.558 (0.015) 13.120 (0.030) 

19 ■■• 16.105 (0.024) 15.412 (0.016) 15.034 (0.022) 14.670 (0.023) 

20 • • • 17.866 (0.031) 17.176 (0.021) 16.781 (0.024) 16.361 (0.025) 



Table 5. Optical Photometry of SN 2000bh 
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^Julian Date minus 2,451,000. 

^1 = ESO NTT 3.6-m. 2 = CTIO 0.9-m. 
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Table 6. Infrared Photometry of SN 2000bh 
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^Julian Date minus 2,451,000. 
^1 = LCO 1-m. 2 = LCO 2.5-m. 



Table 7. Optical Photometry of SN 2000ca 
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^Julian Date minus 2,451,000. 

^1 = ESO NTT 3.6-m. 2 = CTIO 0.9-m. 3 = LCO 1-m. 
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Table 8. Infrared Photometry of SN 2000ca 

JD^' J, H lelescope^ 
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''Julian Date minus 2,451,000. 
H = LCO 1-m. 2 = LCO 2.5-m. 
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Table 9. Optical Photometry of SN 2001ba 
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Table 9 — Continued 

2080.51 ••• ••• 18.261 (0.075) 

2082.49 19.626 (0.078) 18.846 (0.065) 18.339 (0.086) 
2083.51 ••• 18.810(0.114) 

2086.48 • • • 18.698 (0.068) 

2087.49 • • • • • • 18.317 (0.077) 

2088.48 ■ ■ ■ • • • 18.534 (0.109) 
2090.46 20.140 (0.222) 

2091.48 • • • 18.948 (0.098) 18.534 (0.086) 
^Julian Date minus 2,450,000. 



Table 10. Infrared Photometry of SN 2001ba 



JD^ Js H Ks Telescope^ 
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^Julian Date minus 2,450,000. 
^1 = LCO 1-m. 2 = LCO 2.5-m. 



Table 11. K-corrections for Type la Supernovae in LCO Near Infrared Bands^ 



t^\z 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 0.045 0.050 0.060 0.070 0.080 0.090 0.100 0.150 
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-0.122 


-0.146 


-0.171 


-0.194 


-0.239 


-0.286 


-0.339 


-0.407 


-0.485 


-0.878 


19.42 


-0.019 


-0.039 


-0.061 


-0.084 


-0.109 


-0.133 


-0.158 


-0.184 


-0.213 


-0.241 


-0.300 


-0.366 


-0.442 


-0.533 


-0.625 


-1.058 


22.40 


-0.023 


-0.047 


-0.072 


-0.099 


-0.127 


-0.155 


-0.184 


-0.216 


-0.249 


-0.283 


-0.354 


-0.434 


-0.523 


-0.623 


-0.723 


-1.164 


27.46 


-0.034 


-0.069 


-0.104 


-0.137 


-0.173 


-0.206 


-0.242 


-0.280 


-0.321 


-0.361 


-0.448 


-0.540 


-0.642 


-0.752 


-0.862 


-1.319 


31.43 


-0.028 


-0.058 


-0.089 


-0.121 


-0.153 


-0.187 


-0.222 


-0.260 


-0.300 


-0.343 


-0.436 


-0.540 


-0.648 


-0.757 


-0.871 


-1.328 


41.42 


-0.028 


-0.064 


-0.105 


-0.149 


-0.196 


-0.244 


-0.294 


-0.347 


-0.403 


-0.459 


-0.561 


-0.656 


-0.756 


-0.861 


-0.973 




-8.57 


-0.008 


-0.016 


-0.025 


-0.033 


-0.042 


-0.051 


-0.060 


H: 

-0.070 


-0.079 


-0.088 


-0.106 


-0.125 


-0.144 


-0.166 


-0.191 


-0.373 


0.39 


-0.004 


-0.007 


-0.010 


0.007 


0.007 


0.016 


0.004 


-0.001 


-0.006 


-0.013 


-0.030 


-0.049 


-0.077 


-0.108 


-0.146 


-0.397 


1.40 


0.002 


0.005 


0.008 


0.011 


0.013 


0.013 


0.011 


0.007 


0.002 


-0.005 


-0.021 


-0.041 


-0.069 


-0.102 


-0.141 


-0.400 


4.51 


0.008 


0.017 


0.028 


0.038 


0.048 


0.054 


0.057 


0.061 


0.063 


0.065 


0.061 


0.049 


0.024 


-0.007 


-0.045 


-0.290 


8.38 


0.004 


0.011 


0.021 


0.034 


0.047 


0.060 


0.073 


0.087 


0.102 


0.115 


0.133 


0.136 


0.122 


0.098 


0.072 


-0.071 


15.39 


0.015 


0.034 


0.056 


0.082 


0.109 


0.139 


0.170 


0.203 


0.234 


0.263 


0.315 


0.354 


0.384 


0.404 


0.420 


0.549 


19.42 


0.015 


0.035 


0.062 


0.093 


0.126 


0.158 


0.189 


0.218 


0.244 


0.267 


0.306 


0.338 


0.367 


0.393 


0.420 


0.622 


22.40 


0.017 


0.039 


0.065 


0.096 


0.129 


0.162 


0.194 


0.222 


0.249 


0.273 


0.312 


0.346 


0.376 


0.406 


0.436 


0.652 


27.46 


0.011 


0.026 


0.045 


0.069 


0.098 


0.122 


0.144 


0.163 


0.180 


0.196 


0.225 


0.257 


0.281 


0.304 


0.338 


0.607 


31.43 


-0.002 


0.002 


0.012 


0.024 


0.038 


0.051 


0.061 


0.068 


0.075 


0.080 


0.089 


0.097 


0.104 


0.117 


0.145 


0.344 


41.42 


-0.014 


-0.023 


-0.026 


-0.025 


-0.022 


-0.018 


-0.015 


-0.012 


-0.010 


-0.006 


0.001 


0.010 


0.023 


0.046 


0.082 


0.338 


-8.57 


-0.012 


-0.024 


-0.036 


-0.048 


-0.058 


-0.070 


-0.082 


Ks-. 
-0.095 


-0.106 


-0.117 


-0.143 


-0.170 


-0.195 


-0.221 


-0.246 


-0.367 


0.39 


-0.021 


-0.047 


-0.074 


-0.102 


-0.130 


-0.158 


-0.189 


-0.219 


-0.248 


-0.275 


-0.330 


-0.382 


-0.433 


-0.475 


-0.517 


-0.706 


1.40 


-0.025 


-0.053 


-0.082 


-0.112 


-0.143 


-0.175 


-0.208 


-0.239 


-0.269 


-0.300 


-0.362 


-0.421 


-0.477 


-0.530 


-0.580 


-0.813 


4.51 


-0.033 


-0.072 


-0.113 


-0.155 


-0.194 


-0.231 


-0.269 


-0.308 


-0.344 


-0.374 


-0.435 


-0.492 


-0.546 


-0.601 


-0.649 


-0.851 


8.38 


-0.036 


-0.079 


-0.122 


-0.162 


-0.202 


-0.240 


-0.279 


-0.315 


-0.350 


-0.387 


-0.450 


-0.504 


-0.553 


-0.593 


-0.629 


-0.797 


15.39 


-0.0.33 


-0.066 


-0.097 


-0.127 


-0.157 


-0.187 


-0.216 


-0.241 


-0.265 


-0.287 


-0.329 


-0.365 


-0..398 


-0.4.33 


-0.474 


-0.696 


22.40 


-0.018 


-0.035 


-0.051 


-0.066 


-0.080 


-0.092 


-0.099 


-0.105 


-0.106 


-0.108 


-0.117 


-0.130 


-0.148 


-0.179 


-0.216 


-0.534 


27.46 


-0.002 


-0.006 


-0.015 


-0.021 


-0.027 


-0.024 


-0.011 


-0.003 


0.004 


0.008 


0.026 


0.016 


-0.015 


-0.058 


-0.099 


-0.460 


31.43 


-0.004 


-0.012 


-0.024 


-0.036 


-0.043 


-0.043 


-0.036 


-0.025 


-0.015 


-0.006 


0.005 


-0.013 


-0.043 


-0.071 


-0.100 


-0.581 


41.42 


-0.004 


-0.013 


-0.024 


-0.036 


-0.048 


-0.047 


-0.036 


-0.023 


-0.012 


-0.006 


0.007 


0.000 


-0.020 


-0.028 


-0.036 


-0.373 



''The following corrections are subtracted from the photometric data to correct them to the observer's frame (redshift z = 0.00). 
''Time t is the number of days since B-band maximum. 



Table 12. Fitting Supernova Light Curves at Maximum'^ 
Coeff AJ Ai/ AK 



ao 0.080 

ai 0.5104699E-01 

as 0.7064257E-02 

ag -0.2579062E-03 



0.050 
0.2509234E-01 
0.1852107E-02 
-0.3557824E-03 



0.042 
0.2728437E-01 
0.3194500E-02 
-0.4139377E-03 



'^The differential magnitudes with respect to maximum are of 
the form: A mag = ao + S(ai x f), where t is the number of days 
since 5-band maximum, scaled according to the mean of the B- 
and \^-band inverse stretch factors from Jha (2000, Fig. 3.8) 
and divided by {1 + z). The fits are vahd from —12 to +10 d of 
"stretched time" . 



Table 13. V minus Infrared Loci for Slowly Declining Type la Supernovae'' 



Paranieter \ color index 


V - J 


V - H 




V - K 




ao 


-0.902(0.005) 


-1.194 




-0.971 




ai 


-0.07136(0.00104) 


-0.4793682E- 


-01 


-0.4089032E- 


-01 






0.4647584E- 


-02 


0.3613656E- 


-02 


as 




0.8813858E- 


-04 


0.6178713E- 


-04 


a4 




-0.5205402E- 


-05 


-0.3648954E- 


-05 


Range^ 


[-8,+9.5] 


[-8,+8.5] 




[-8,+27] 




RMS residual 


± 0.067 


± 0.062 




± 0.138 






3.5 


1.9 




4.9 





^The range of Amis (5) is 0.81 to 1.00. The fits are of the form: V - [J, H, K] = ao 
+ E(aj X f), where t is the number of days since S-band maximum, corrected for time 
dilation but not stretch. 

''The range in days with respect to B-band maximum for which the fit is most valid. 
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Table 14. V minus Infrared Loci for Type la Supernovae of Mid-Range Decline Rates^ 



Parameter \ color index 


V-H 




V-K 




ao 


-0.951 




-0.741 




ai 


-0.5383673E- 


-01 


-0.5071254E- 


-01 


a2 


0.5913753E- 


-02 


0.3281215E- 


-02 


as 


0.1155147E- 


-03 


0.2158419E- 


-03 


04 


-0.7144163E- 


-05 


-0.8056540E- 


-05 


Range^ 


[-9,+27] 




[-9,+27] 




RMS residual 


± 0.093 




± 0.105 






4.2 




3.9 





'^Krisciunas et al. (2000, Tabic 9) gave hinged two-line fits to 
the color curves. Actual data are redder at t = +6 d than the two- 
line fits. Here we give fourth order polynomial fits to the same 
extinction-corrected data. The fits are of the form: V — [J, H, K] 
— ao + T:{ai X f), where t is the number of days since S-band 
maximum. 

^The range in days with respect to S-band maximum for which 
the fit is most valid. 
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Fig. 1. — Near infrared photometry of SN 1999ee. For J, and H the Las Campanas data are 
red filled in symbols, while the YALO data are filled in with blue. The YALO data include 
filter-based corrections to place them on the system of Persson et al. (1998) used at Las 
Campanas Observatory. 

Fig. 2. — Finding chart for SN 2000bh. North is up and east to the left. The image is 6.5 
arcmin on a side. It is a V-band image obtained on the CTIO 0.9-m telescope. The SN is 
indicated by the short line. 

Fig. 3. — Optical and infrared photometry of SN 2000bh. We have added BVI fits derived 
using the Amis method of Phillips et al. (1999). 

Fig. 4.— Similar to Fig. 2, but for SN 2000ca. 

Fig. 5. — Optical and infrared photometry of SN 2000ca. 

Fig. 6. — Similar to Fig. 2, but for SN 200 Iba. This V-band image was obtained with the 
Las Campanas 2.5-m telescope. 

Fig. 7. — Optical and infrared photometry of SN 2001ba. 

Fig. 8. — J-band light curves of 8 Type la SNe, normalized to maximum brightness and 
stretched in the time axis according to inverse stretch factors derived from Fig. 3.8 of Jha 
(2002). We have also take out the time dilation. The key to the symbols is: 1980N (solid 
dots), 1986G (large open squares), 1998bu (larger open circles), 1999aw (X's), 1999ee (LCO 
data = small open squares, corrected YALO data = triangles), 2000ca (pinched in squares), 
2001ba (stars), 2001el (smaller open circles). 

Fig. 9. — JH maxima of 8 Type la SNe and X-band maxima of 6 Type la SNe. The time 
axis is "stretched days" (see Fig. 8 and the text). We also give third order polynomial fits. 

Fig. 10.— /-band data of SNe 1998bu, 1999aw, 1999ee, 2000ca, 2001ba, and 2001el. The 
data are stretched in the time axis in a manner similar to Figs. 8 and 9. 

Fig. 11.— V minus IR colors of Type la SNe with dechne rates 0.81 < Ami5(5) < 1.00. The 
data are K-corrected, and have been corrected for extinction due to dust and time dilation, 
but not stretch. The solid line in the top plot shows the range of uniformity oiV — J . The 
dashed lines in the bottom two plots are the loci for mid-range decliners (Krisciunas et al. 
2000, 2003). The solid lines in the bottom two plots are fourth order polynomial fits to the 
data from —8 to -\-21 days with respect to S-band maximum. On the whole the random 
errors of the photometry are on the order of the size of the points, or smaller. In the cases 
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of SNe 1999gp and 1999aw we show the considerably larger error bars. The points are color 
coded as follows: SN 1999aa (magenta); 1999aw (red), 1999ee (blue), 1999gp (cyan), 2000ca 
(green), and 2001ba (yellow). 

Fig. 12. — Dereddened colors of Type la SNe vs. the dechne rate parameter Ami5(i?). The 
(green) triangles correspond to objects which have few infrared observations in the [—12, +10] 
day window, so their IR maxima are not as well determined; these objects were not used for 
the regression lines shown. 

Fig. 13. — Dereddened V — J and V — H colors of Type la SNe at 6 days after the time of 
the S-band maxima, vs. the decline rate parameter Amis (5). 

Fig. 14. — R — z vs. V — I colors for spectrophotometric standards (Stone & Baldwin 
1983; open circles) and for Landolt (1992a) standards. The z magnitudes of Stone-Baldwin 
standards were taken from Hamuy (2001, Appendix B); the solid line is a fit to data for 
these stars. The dashed hue is a fit to points with V — I > 0.97. Three particular stars are 
Rubin 149C (open squares), SA95-96 (diamond), and Rubin 152E (X's). See text for further 
details. 
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Krisciunas et al. Fig. 2 
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Krisciunas et al. Fig. 6 
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Time since B— band maxinnunn (stretched days) 



Krisciunas et al. Fig. 8 
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Time since B-band maximum (stretched days) 



Krisciunas et al. Fig. 10 
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Krisciunas et al. Fig. 11 
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Krisciunas et al. Fig. 14 



